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Key Objectives 


¢ Continuation of previous work, which: 
— Compared flashes generated by flash rate parameterization schemes 
(FRPSs) ina WRF-Chem model simulation with lightning observations: 
¢ Oklahoma Lightning Mapping Array (OK LMA) 
¢ National Lightning Detection Network (NLDN) 
— Tentatively concluded lightning-generated NO, (LNO,) production is 
around 125 moles flash* 


¢ Current work objectives: 

— Define and incorporate new lightning flash channel vertical 
distributions and IC:CG ratios into the WRF-Chem model based on 
lightning data from a LMA for the storm of interest 

— Analyze distribution of observed and model-simulated trace gas 
species in storm inflow and outflow 

— Determine NO production scenario for IC and CG lightning 


Severe convection developed ~21Z May 


Background 


29 along KS/OK border and continued 
until 04Z May 30 


Aircraft sampled storm and its 
environment from 20Z May 29 to 01Z 
May 30 


DC-8 focused on storm inflow & outflow 
GV & Falcon concentrated on outflow 


Ground-based data included: 


Dual-Doppler radar (NEXRAD level II 
regional) 


Shared Mobile Atmospheric Research 
and Teaching Radar (SMART-Radar) 


NLDN cloud-to-ground flash data 
OK LMA flash initiation density data 





NEXRAD Composite Reflectivity 2240Z on 29 May 
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Blue circles: LMA stations 
Green outline: Extent of 3-D lightning mapping capability 
Gray outline: Extent of 2-D lightning detection 3 


WRF-Chem Model V3.6.1 


¢ Grid resolution: dx = dy = 1-km, dz = 50-250 m 
¢ Initialized with 18Z NAM ANL (6-hr) for boundary conditions 
¢ Lightning Data Assimilation from 18-212Z (Fierro et al., 2012) 


Microphysics 

Planetary boundary layer 
Land surface 

Radiation (short & longwave ) 
Photolysis 

Trace gas chemistry 


Flash rate 


LNO 


Morrison 

Yonsei University (YSU) 

Noah 

Rapid radiative transfer model for GCMs (RRTMG) 
F-TUV 

MOZART 


> Updraft volume based on AL supercells (UP510_S) 
> Coarsely prescribed IC:CG ratios (Boccippio et al., 2001) 
replaced with IC:CG ratios based on LMA and NLDN obs 


Flash segment vertical distribution based on observations 


LNO, Parameterization Scheme 


Replaced the typical lightning flash channel 
distributions (DeCaria et al., 2000; 2005) with 
observed IC & CG vertical distributions 
— Used flash channel segment data from observed 
storm’s respective LMA network 


— IC &CG distributions for 29 May both appear to 


be single Gaussian where channels maximize at New IC & CG 
~10km (-42°C) distributions for 

— Previous distributions were set to maximize the Mee) Bs 29 May case 
lightning channels at -15°C (IC & CG) and -45°C © 2) (OS km") F(z) (0.5 km") 


(IC), or 6km and 10.5km, respectively 








Found 125 moles flash? provided best fit with 
observed anvil NO, when using DeCaria et al. 
vertical distributions. This scenario: 


—  |Is much smaller than mean value of 500 moles 
flash? found in previous mid-latitude simulations 
(Ott et al., 2010) 


— Will be tested using new distributions 





Horizontal placement of NO based on reflectivity 
> 20 dBZ in each grid cell 


DeCaria et al. 
(2000, 2005) 











0.0 


0.2 
f(z) (km*') 








Max vertical velocity (Wmax) 
Cloud top height (CTH) 


Updraft volume (UpVo!) 


Ice water path (/WP) 
Precipitation ice mass (P/M) 

Ice mass flux product (/MFP) 
Graupel volume (CSU_ GEV) 
35-dBZ volume (CSU_VOL35) 
Precipitation ice mass (CSU_ PIM) 


Graupel echo volume 
(-40°<T<-5°C; ALGEV5) 


ALGEV5 for supercells (ALGEV5_S) 


Graupel echo volume 
(-40°<T<-10°C; ALGEV10) 


ALGEV10 for supercells 
(ALGEV10_S) 


Updraft volume(w>5m s+, 
-40°C<T<-10°C; ALUP510) 


ALUP510 for supercells 
(ALUP510_S) 


Methodology 


¢ Compared flash rate trends over the 
observed and model-simulated storm’s 
lifetime 


| Reference 


Price & Rind (1992) 
Price & Rind (1992) 


Deierling & Petersen 
(2008) 


Petersen et al. (2005) 
Deierling et al. (2008) 
Deierling et al. (2008) 
Basarab et al. (2015) 
Basarab et al. (2015) 
Basarab et al. (2015) 


Carey et al. (2015) 


L. Carey 


Carey et al. (2015) 


L. Carey 


Carey et al. (2015) 


L. Carey 


Used 15 different FRPS, including those 
from literature, as well as recently 
developed schemes from DC3 radar and 
LMA data 


Selected the FRPS that reasonably 
represented the total observed flashes 
over the storm and the flash rate trends 
from the LMA 


¢ Assumed LNO, production is 125 moles 
flash* and will adjust as necessary 


¢ Analyzed trace gas species (i.e., CO, NO,, 
O,) using model-simulated values and 
aircraft (DC-8 & GV) observations to: 


Create probability distribution function 
(PDF) plots in storm outflow 


Evaluate convective transport 


Determine best fit NO production 
scenario 


IC:CG Ratios 


NLDN CG & LMA TOTAL FLASHES 21:10—04:192Z 





LMA Total Flashes: 31633 


¢ Coarsely prescribed IC:CG ratios sooof. NLON CG Flashes: 8679 
from Boccippio et al. (2001) 
provide a mean IC:CG of 3.90 + 
0.49 over the region where the 
severe convection occurred 


¢ LMA total and NLDN CG flashes 
indicate the IC:CG ratio fluctuates le 
over the lifetime of the i a 
convection on 29 May 10—MINUTE IC:CG RATIOS 21:10—04:19Z 


— Mean I!1C:CG ratio over storm 
lifetime is 2.73 +2.51 


2000 


FLASHES 10—MINUTE™ 


1000 


v 4 > 


¢ Time evolving IC:CG ratios are 
applied in the model to the storm 
of interest, while climatological 
IC:CG values are used in the area 
Surrounding the storm 


IC:CG RATIO 10—MINUTE™ 





Model Flash Rates vs. Observations 


MODEL STORM START TIME ADJUSTED 40MIN EARLY 


MODEL STORM START TIME ADJUSTED 40MIN EARLY 
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Time (UTC) 
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Total Flashes 10—Minutes™ 
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496436 
471270 
18062 
393313 
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MP 97: 95:0 ga 950° 


Observed storm occurs 21:10-04:10 UTC 


Model-simulated storm onset is delayed 
40 min (21:50 UTC) 


Plotted observed & model-simulated 
flash rates to both begin at 21:10 UTC 


Model Flash Rates vs. Observations 


¢ Selected the updraft volume FRPS 
for Alabama supercells (ALUP510_S) 
for use in model 
— Based on updraft (w>5 ms") volume 


and mixed-phase region (-40°C <T < 
-10°C) 

r UpVol (—): 32340 
00 F- LMA OBS (—): 18911 


¢ All FRPSs incorporating 
hydrometeors overestimate (by 
7-74x) the total flash observations 


¢ All FRPSs developed for the Colorado 
region overestimate (by 7-11x) the 
total flash observations 
— 35-dBZ volume was generally 
developed from storms with shallow 
warm cloud depths (< 1km), while the 


29 May case has a deep warm cloud 
depth (2.5km) 
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FREQUENCY OF OBSERVATIONS 


CO & O; PDFs in Storm Outtlow 


OUTFLOW CO PDF BETWEEN 10.23-10.48 km ALTITUDE OUTFLOW CO PDF BETWEEN 10.48-10.72 km ALTITUDE 
1.0 1.0 
0.9 LEGEND: 0.9 LEGEND: 
mepcs ( 244) mcy (:«‘177) 
0.8 MB MODEL (21698) 0.8 HE MODEL (12639) 





Model-simulated 
CO peaks at larger 


0.7 


0.6 


0.5 


FREQUENCY OF OBSERVATIONS 





. 0. values than 

0.5 0.3 observations 
0.2 0.2 

0.1 0.1 

0.0 0.0 

:,rft Pr FF FCF Sf LY fF LH e 
CO (ppbv) CO (ppbv) 
OUTFLOW QO, PDF BETWEEN 10.23-10.48 km ALTITUDE | OUTFLOW QO, POF BETWEEN 10.48-10.72 km ALTITUDE 


LEGEND: 
mepcs ( 244) 
@= MODEL (21698) 


LEGEND: 
me GV ( 177) 
MB MODEL (12639) 


Model-simulated 
O, peaks at 
smaller values 
than observations 


FREQUENCY OF OBSERVATIONS 
FREQUENCY OF OBSERVATIONS 





© 9 5 © © 
© : \’ xb xv 


wePe Cf Ff pg oe ££ 
O; (ppbv) O; (ppbv) 


* Comparison of peak vertical velocities (not shown) suggests model-simulated values may be 
20% greater than the SMART-Radar observations 10 


NO. PDEs in Storm Outflow 


OUTFLOW NO, PDF BETWEEN 10.23-10.48 km ALTITUDE OUTFLOW NO, PDF BETWEEN 10.48-—10.72 km ALTITUDE OUTFLOW NO, PDF BETWEEN 10.72-—10.97 km ALTITUDE 
‘ 1.0 1.0 





LEGEND: 0.9 LEGEND: 0.9 LEGEND: 
mDc8 ( 103) MGV 6 (177) mm GV ( 157) 
HB MODEL (21973) 0.8 @B MODEL (12474) 0.8 MB MODEL (12504) 
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FREQUENCY OF OBSERVATIONS 

FREQUENCY OF OBSERVATIONS 
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Conclusions 


A single model domain at fine resolution (1-km) produces a storm 
of roughly the same size and with similar characteristics as the 
observed 


FRPSs based on hydrometeors or Colorado storms are not ideal 
for the severe Oklahoma convection observed on 29 May 


— Selected a FRPS based on updraft volume (w > 5 ms“) within the 
mixed-phase region of Northern Alabama supercells 


LNO, production may be closer to 250 moles flash* 


— Removed influence of an IC upper lightning channel on NO, by 
replacing the IC & CG vertical distributions (DeCaria et al. 2000, 2005) 
with observed lightning channel distributions from 29 May 


— May be less than the mean value for mid-latitude storms (500 moles 
flash) due to presence of smaller flashes 


Future Work 


29-30 May 2012 Oklahoma severe convection: 
— Finalize NO production scenario for IC and CG LNO, scheme 


— Test other NO production scenarios from the literature and compare against the 
scenario selected for 29 May 
¢ 500 moles flash? (Ott et al., 2010) 
¢ Lightning Nitrogen Oxides Model (LNOM) results for IC & CG flashes (Koshak, 2014) 


— Investigate O; changes within the cloud and downwind of the storm 


6-7 June 2012 Colorado squall line: 
— Test and select one of the 15 FRPSs in the WRF model 
— Determine NO production scenario for IC and CG LNO, scheme 
— Investigate O, changes within the cloud and downwind of the storm 


Compare the results between the storms to: 
— Investigate which FRPSs are most appropriate for the two types of convection 
— Examine the variation in LNO, production 


Compare simulated LNO, results from WRF-Chem with other previously 
studied mid-latitude thunderstorms 
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